Dynamics of Brane- World Cosmological Models 

A. A. Coleyi'* 

^Department of Mathematics and Statistics, Dalhousie University, Halifax, Nova Scotia 

We show that genericaUy the initial singularity is isotropic in spatially homogeneous cosmological 
models in the brane-world scenario. We then argue that it is plausible that the initial singularity 
is isotropic in typical brane world cosmological models. Therefore, brane cosmology naturally gives 
rise to a set of initial data that provide the conditions for inflation to subsequently take place, 
thereby solving the initial conditions problem and leading to a self-consistent and viable cosmology. 

98.80.Cq 



O 
O 

(N 

Oh; 

< 

in 

> 

^ , 
o . 
o ■ 



o : 

Oh! 



X 



I. INTRODUCTION 

Recent developments in string theory suggest that 
gravity may be a truly higher-dimensional theory, becom- 
ing effectively 4-dimensional at lower energies. This leads 
to modifications to Einstein's theory of general relativ- 
ity (GR) at high energies and particularly at early times. 
There is currently great interest in higher-dimensional 
gravity theories inspired by string theory in which the 
matter fields are confined to a 3-dimensional 'brane- 
world' embedded in 1 -I- 3 -I- c? dimensions, while the grav- 
itational field can also propagate in the d extra dimen- 
sions (i.e., in the 'bulk') In this paradigm it is not 
necessary for the d extra dimensions to be small or even 
compact, a departure from the standard Kaluza-Klein 
scenario. In recent work Randall and Sundrum |^ have 
shown that for d = 1, gravity can be localized on a sin- 
gle 3-brane even when the fifth dimension is infinite. An 
elegant geometric formulation and generalization of the 
class of Randall-Sundrum-type brane-world models has 
been given in 

The dynamical equations on the 3-brane differ from 
the GR equations by terms that carry the effects of 
imbedding and of the free gravitational field in the 5- 
dimensional bulk. The local (quadratic) matter fields 
corrections are significant only at very high energies. In 
addition, there are nonlocal effects from the free gravi- 
tational field in the bulk, transmitted via the projection 
£f^i, of the bulk Weyl tensor, that contribute further cor- 
rections to the Einstein equations on the brane. can 
be irreducibly decomposed (with respect to a timelike 
congruence u^^) in terms of an effective nonlocal energy 
density on the brane, U, arising from the free gravita- 
tional field in the bulk, an effective nonlocal anisotropic 
stress, T'p,i,, on the brane, and an effective nonlocal en- 
ergy flux on the brane, [Q. 

In general, the conservation equations do not deter- 
mine all of the independent components of S^^i, on the 
brane; in particular, there is no evolution equation for 
Vfj^i,. Thus in general, the projection of the 5-dimensional 
field equations onto the brane does not lead to a closed 
system. If the induced metric on the brane is flat, and the 



bulk is anti-de Sitter, as in the original Randall-Sundrum 
scenario then = 0. More importantly, in cos- 
mology the background induced metric is not flat, but 
a spatially homogeneous and isotropic Robertson- Walker 
(RW) 'Friedmann' model, for which 



0, 



(1) 



where is the totally projected part of the brane covari- 
ant derivative. Hence, in such cosmological settings the 
evolution of is fully determined, and the system of 
equations is closed. In general U ^ in the Friedmann 
background. Much effort is currently being devoted to 
understand the cosmology of the brane world scenario, 
and the Friedmann brane models have been extensively 
investigated |^,D). 

Cosmological observations indicate that we live in a 
Universe which is remarkably uniform on very large 
scales. However, the spatial homogeneity and isotropy 
of the Universe is difficult to explain within the stan- 
dard GR framework since, in the presence of matter, 
the class of solutions to the Einstein equations which 
evolve towards a RW universe is essentially a set of mea- 
sure zero [0. In the inflationary scenario, we live in 
an isotropic region of a potentially highly irregular uni- 
verse as the result of an expansion phase in the early 
universe thereby solving many of the problems of cos- 
mology. Thus this scenario can successfully generate a 
homogeneous and isotropic RW-like universe from initial 
conditions which, in the absence of inflation, would have 
resulted in a universe far removed from the one we live 
in today. However, still only a restricted set of initial 
data will lead to smooth enough conditions for the onset 
of inflation (i.e., the so-called cosmic no-hair theorems 
only apply to non-generic models Q ) , so the issue of ho- 
mogenization and isotropization is still not satisfactorily 
solved. Indeed, the initial conditions problem, that is 
to explain why the Universe so isotropic and spatial ho- 
mogeneous from generic initial conditions, is one of the 
central problems of modern theoretical cosmology. We 
would like to revisit these issues in the context of brane 
cosmology. 
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A. Governing equations 

The field equations induced on the brane, using the 
Gauss-Codazzi equations, matching conditions and Z2 
symmetry, were derived in [^[^], resulting in a modifica- 
tion of the standard GR equations with the new terms 
carrying bulk eff'ects onto the brane: 

G = ^9^V + + K l5'p;y £ = hi T , (2) 

where — Stt/M^jX = 6k^/k'^. The generalized Fried- 
mann equation, which determines the expansion of the 
universe, in the case of spatially homogeneous cosmolog- 
ical models is given by 

9 In/ p \ lo It 1. 2U 

where is the scalar curvature of the hypersurfaces or- 
thogonal to the fluid flow, which we associate with the 
cosmological fluid, 2a^ = (T°'''(Tab is the shear scalar, and 
H = a/ a is the mean Hubble expansion parameter. 

The brane energy-momentum tensor for a perfect fluid 
or a minimally- coupled scalar field is given by 

Tf^u = pu^_,Uu + phf^„ . (4) 

where is the 4- velocity, p and p are the energy density 
and isotropic pressure, and h^i, = g^u -\- Uf^u^. It follows 
1^ that the brane energy-momentum tensor separately 
satisfies the conservation equation, 

p + 3H{p + p) = 0. (5) 

For a minimally coupled scalar field the energy density 
and pressure are, respectively, 

P=^4>' + V{c^), p=^<P^-V{^), (6) 

and Eqn (^) is equivalent to the Klein-Gordon equation. 

There are many reasons to consider cosmological mod- 
els that are more general than RW, both spatially homo- 
geneous and anisotropic, and spatially inhomogcncous. 

(i) The 3-curvature in RW models is given by = 
6k/a^, where k = 0, ±1 is the curvature constant. An 
equivalent 3-curvature occurs in spatially homogeneous 
and isotropic curvature models, and a similar term occurs 
in other cosmological models. 

(ii) For BRW models, equation (|l|) implies that lA — 
U{t). In the models of Bianchi type I, Qa = but there 
is no restriction on Vab] however, the particular case in 
which this term is zero can be studied. Thus, in RW and 
Bianchi I models the evolution equation for W is Q 

U + AHU = , (7) 

which integrates to U — Uo/a'^, which has the structure 
of a 'dark' radiation fluid, where Uq can be negative. 



(iii) A Bianchi I brane is covariantly characterized in 
1^. The conservation equations reduce to Eqn (||), an 
evolution equation for U and a differential constraint on 
■ppiy. The presence of in the governing equations 
means that in general we cannot integrate to find the 
shear as in GR. However, when the nonlocal energy den- 
sity vanishes or is negligible, i.e., U ~ 0, then the conser- 
vation equations imply cr'^'^V^n, — 0, which is consistent 
on the brane Q. This assumption is often made in the 
case of RW branes Q . The shear evolution equation may 
then be integrated to give 

't'^^'t,.. = — , S = 0. (8) 

Bianchi I models on the brane have been studied by a 
number of authors |^,^ . A similar shear term occurs in 
other hypersurface orthogonal Bianchi cosmological mod- 
els (such as Bianchi type V models). 

(iv) We assume that the matter content is a non-tilting 
perfect fluid with a linear barotropic equation of state, 
i.e., P = (7 — where the energy conditions imply 
p > 0, and the constant 7 satisfies 7 G [0, 2]. A dynami- 
cal analysis of scalar field models indicates that at early 
times the scalar field is effectively massless. A massless 
scalar field is equivalent to a perfect fluid with a stiff 
equation of state parameter 7 = 2. Eqn (||) then yields 
p — poa~^"', where pQ > 0. 

We can therefore write down a phenomenological gen- 
eralized Friedmann equation: 

^ -3^ + 6A^"^+3'^+^ + ^- 

This equation is applicable in a wide class of spatially ho- 
mogeneous cosmological models; in particular, it is valid 
in Bianchi I brane models. In many applications the 
four-dimensional cosmological constant is assumed to be 
zero 1^; here we shall assume that if it is non-zero it is 
positive, i.e. A > 0. 

II. EARLY UNIVERSE: 

It is of considerable interest to study the classical dy- 
namical effects in these cosmological models. A unique 
feature of brane cosmology is that p^ dominates at early 
times which will lead to completely different behaviour 
to that in GR. First, there is the question of the ex- 
istence of singularities. The generalized Raychaudhuri 
equation governs gravitational collapse and initial singu- 
larity behaviour on the brane. The local energy density 
and pressure corrections further enhance the tendency to 
collapse (if 2p + 3p > 0) jlO|. The nonlocal term can 
act either way depending on its sign; a negative U en- 
hances the localization of the gravitational field on the 
brane (the singularity can be avoided in this case), and 
a positive U acts against localization, and also reinforces 
the tendency to collapse. For the models governed by 



2 



Eqn it is easy to show that once a critical value of a 
is attained, a singularity must occur. 

As in GR the powerful singularity theorems of Penrose 
and Hawking, which guarantee the existence of spacetime 
singularities, will be generally applicable in brane world 
models. However, these theorems give little information 
about the nature of the singularities they predict. The 
most detailed proposal for the structure of spacetime sin- 
gularities in GR are the conjectures of Belinskii, Khalat- 
nikov and Lifshitz (BKL) which essentially consist of 
two parts. The first is: I. Perfect fluid spacetimes (with 
a linear equation of state) with non-stiff matter have 
the property that asymptotically close to the singularity, 
matter is not dynamically significant. In the case of spa- 
tially homogeneous models this implies that spacetimes 
are space-like and oscillatory (asymptotic Mixmaster be- 
haviour, oscillating indefinitely as the cosmological initial 
singularity is approached into the past), or are space- 
like and non-oscillatory (and asymptotically Kasner at 
the singularity). In the case of stiff matter, which in- 
cludes the massless scalar field case, the matter is not in- 
significant near the singularity and generically the space- 
times have singularities which are space-like and non- 
oscillatory (asymptotically of Jacobs form). 

Theoretical justification for part I of the BKL con- 
jectures in Bianchi models, in which the Einstein equa- 
tions are a system of ODE's. has been provided using dy- 
namical systems methods [ [l2| and more rigorous meth- 
ods [ p3[ |. There are, however, special classes of models 
that do not obey the BKL conjectures in GR, the most 
important being models with an isotropic initial singu- 
larity fl^ , whose evolution near the cosmological initial 
singularity is approximated (in an appropriately defined 
mathematical sense) by the flat Friedmann model. 

An essential feature of the brane models governed by 
Eqn (^) is that at high densities the term {n'^ p^/GXa^^) 
dominates and the the effective equation of state becomes 
ultra stiff. Consequently matter dominates the shear and 
curvature (and the other) terms in Eqn (^ when 7 > 1, 
leading to isotropic expansion of the early universe in 
such cases. In particular, a^/H'^ ^ as t — > 0. 

Indeed, the flat, spatially homogeneous and isotropic 
non-GR brane-world (without brane tension) BRW 
model, denoted here by J^b, in which a{t) ~ ||], is 
always a source/repeller for 7 > 1 (as can clearly be 
seen from Eqn (^). These BRW models are valid as the 
initial singularity is approached (t — > 0), and therefore 
for all physically relevant values of 7 the singularity is 
isotropic [0. We expect this to be a generic feature 
of more general cosmological models in the brane-world 
scenario, as we will discuss this further below. 

From Eqn (|^), in the brane-world scenario anisotropy 
dominates at early times only for 7 < 1 (whereas in GR 
it dominates for 7 < 2!), in which case the repellers are 
the usual anisotropic Kasner models. In the absence of 
shear, the BRW models are sources for 7 > 1/3 when 
and 7 > 2/3 when U ^ The fact that 



the initial singularity is isotropic in Bianchi type I and V 
models was noted in [^^. 

The particular models we have considered, in which 
the curvature and shear are given by the expressions in 
the above phenomenological Eqn (|9|), are special. In par- 
ticular, the Bianchi I and V models are not generic, and 
so the study of the dynamics of these models does not 
shed light on the typical behaviour of spatially homoge- 
neous brane models. In order to do this we shall next 
consider the general Bianchi type IX (paradigm) model. 

The matter corrections to the Einstein equations on 
the brane are given by 



•Spi/ = j2P^'^f^^'^ + T2/0 (P + 2p) h 



for a perfect fluid (or minimally-coupled scalar field). 
There are also nonlocal effects from the free gravitational 
field in the bulk via lA. All of the bulk corrections may 
be consolidated into an effective total energy density via 
Eqn (^, so that the modified Einstein equations take on 
the standard perfect fluid form with 



1^6 



Ay 



A 



^ 2A A 3k6 k2 



For the Bianchi IX models with 7^ = diag{a\, a^, flg), 
where the mean scale factor is given in terms of the three 
spatial scales by a^ = aia2az, it can be proven [p^ from 
the conservation Eqn (^ , the nonlocal conservation Eqn 
(0) , and the perfect fluid generalized Raychaudhuri equa- 
tion that there is a singularity (at t — Q), and that p — > 00 
as t ^ 0. Close to singularity we have that H < 0, and 
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The Bianchi IX equations of motion then become 



^^ 2\ll / 2 2\2 

(InaJ) — (a2 — flq ' 



^3 7 



«^(l-7) 
A 



a^p^, et eye. (11) 



where a prime denotes differentiation with respect to r 
(where dt = a^dr so that r — > —00 as the singularity 
is approached), and et eye. denotes two more equations 
obtained by cycling the indices of Oi on the left-hand-side 
of Eqn (the right-hand-side remains the same). As 
in the GR case, there exists a first integral 



(Ina^)'(lna^ 



2alal-aj + cyc. = —a''p\ (12) 



The BRW solution J-^, ai = 02 = ^ a = Ob (with 
ab{t) = t~), occurs as r — > —00 (t 0) and a — > 00. 
From the governing equations (where we use Eqn ( |T^ ) 
to eliminate K^p^/X in Eqns (|0) and the Raychaudhuri 
equation governs the evolution of H), we obtain the fol- 
lowing 5 essential eigenvalues for the linearization about 
J'b- 

3(7 - 1), 3(7 - 1), (37 - 1), (37 - 1), (37 - 1) (13) 
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(the first two eigenvalues correspond to shear modes 
while the last three correspond to curvature modes). 
Thus !Fb is a past-attractor (repeller/source) of the 
Bianchi type IX models for 7 > 1 (for 7 = 1 the BRW 
solution is an analogue of the Jacobs stiff fluid solution 
in GR and is a past-attractor) . There are no other past- 
attractors in particular, the analogues of the non-flat 
Kasner vacuum solutions are saddles. Hence the singu- 
larity is non-oscillatory and isotropic. It can be shown 
that J-f, is a source in all Bianchi models |0 . 

The results are incomplete in that a description of the 
gravitational field in the bulk is not provided. Unfor- 
tunately, the evolution of the anisotropic stress part is 
not determined on the brane. These nonlocal terms also 
enter into crucial dynamical equations, such as the Ray- 
chaudhuri equation and the shear propagation equation, 
and can lead to important effects. The correction terms 
must be consistently derived from the higher-dimensional 
equations. Additional modifications occur for more gen- 
eral higher-dimensional (bulk) geometry, higher curva- 
ture corrections, higher-dimensional matter fields in the 
bulk, and for motion of the brane. 



III. DISCUSSION 

Consequently, we have shown that generically the ini- 
tial singularity is isotropic in spatially homogeneous 
brane world cosmological models. 

The second part of the BKL conjectures is: II. Each 
spatial point evolves towards the singularity as if it were a 
spatially homogeneous cosmology. That is, generic space- 
times have the property that spatial points decouple near 
the singularity and the Einstein equations effectively re- 
duce to ODEs, so that the local dynamical behaviour is 
asymptotically like Bianchi models near the singularity. 
Therefore, according to the BKL conjectures, the sin- 
gularities in general four-dimensional spacetimes in GR 
are space-like and oscillatory or are space-like and non- 
oscillatory (e.g., for massless scalar fields). 

An analysis of the behaviour of spatially inhomoge- 
neous solutions to Einstein's equations near an initial sin- 
gularity is in its infancy, and hence there is less support 
for part II of the BKL conjectures. However, a special 
class of Abelian G2 spatially inhomogeneous models were 
analysed and it was found that the evolution at differ- 
ent spatial points approach that of different Kasner solu- 
tions . A recent numerical investigation of a class of 
vacuum Gowdy G2 cosmological spacetimes has shown 
evidence that at a generic point in space the evolution 
towards the initial singularity is asymptotically that of 
a spatially homogeneous spacetime with Mixmaster be- 
havior In both of these cases the presence of the 
inhomogeneity ceases to govern the dynamics asymptot- 
ically toward the singularity, thereby providing further 
support for the BKL conjectures. 

In addition, in a recent qualitative analysis of a class 



of spatially inhomogeneous G2 brane cosmological mod- 
els (with one spatial degree of freedom) near the initial 
cosmological singularity, it was found that J-b is again a 
local source It was also argued, based upon local dy- 
namical considerations and physical arguments, that the 
main result that the singularity is isotropic will persist 
when additional, more general, affects are included. 

Thus it is plausible, from the BKL conjectures in 
GR and from this recent study of spatially inhomoge- 
neous brane models, that typically the initial singularity 
is isotropic in brane world cosmological models. 

Therefore, unlike the situation in GR, it is plausible 
that a wide range of brane cosmological models (and of 
non-zero measure) admit an isotropic singularity. Such 
a 'quiescent' cosmology po| , in which the universe be- 
gan in a highly regular state but subsequently evolved 
towards irregularity, might offer an explanation of why 
our Universe might have began its evolution in such a 
smooth manner and may provide a realisation of Pen- 
rose's ideas on gravitational entropy and the second law 
of thermodynamics in cosmology 1 2W . 

More importantly, it is therefore possible that a qui- 
escent cosmological period occuring in brane cosmology 
provides a physical scenario in which the universe starts 
off smooth and that naturally gives rise to the conditions 
for inflation to subsequently take place. As noted ear- 
lier, in the conventional scenario only a restricted set of 
initial data will lead to the conditions possible for infla- 
tion to occur. We thus argue that brane cosmology in 
tandem with inflation may lead to a fully self-consistent 
and physically viable cosmology Q . 

Subsequently, the models evolve essentially as in stan- 
dard cosmology. The modifled Friedmann equation for a 
flat RW brane (with A = = U) and the Klein-Gordon 
Eqn (||) yield p = —34)^p at early times for a scalar field 
source, so that p is monotonically decreasing and the 
models will eventually evolve to the low density regime. 

The intermediate dynamics is affected by the brane 
corrections. The issue of inflation on the brane was inves- 
tigated in [|o| , where it was shown that on an RW brane 
in 5-dimensional anti de Sitter space the quadratic term 
in p increases friction in the inflaton field equation and 
infiation at high energies proceeds at a higher rate than 
the corresponding rate in GR. Moreover, it was shown 
that, contrary to expectations, a large initial anisotropy 
introduces more damping into the scalar field equation 
of motion and results in more infiation |^ . 

At late(r) times a number of features of the cosmo- 
logical dynamics can be deduced directly from Eqn (||), 
consistent with the qualitative analysis of perfect fiuid 
RW and Bianchi type I and V cosmological models in 
the Randall-Sundrum brane- world scenario of H. Mod- 
els with a positive curvature can recoUapse. However, for 
U < models can (re) collapse (even without a positive 
curvature) for any values of 7. Indeed, for U < and 
positive curvature (as in the k = 1 BRW models), there 
exist oscillating universes in which the physical variables 
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oscillate per iodically without reaching any spacelike sin- 
gularity |15|] (it was noted earlier that when there are 
bulk effects present a singularity can be avoided) . When 
a positive cosmological constant is present, the de Sit- 
ter model is always the global attractor for U > 0. For 
U < 0, models can (re)collapse (even without a positive 
curvature) so that in this case the de Sitter models is 
only a local attractor (and the cosmic no-hair theorem is 
consequently violated @). 
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